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Abstract 



CIhI We present an implementation of heavy quarkonium production within 



a perturbative QCD cascade based on the Color Dipole Cascade model. 
We consider the processes most relevant in the context of the ip' sur- 
plus at the Tevatron; g ^ ip' and c ^ m. the color-singlet model 
^ ' and g ^ Tp' through the color-octet mechanism. Our implementation 

■ is, however, easily extendible to other quarkonia and other production 

mechanisms. 

Where comparison is possible we find good agreement with ana- 
lytical calculations. 

We present some suggestions for measurements at the Tevatron 
that would be sensitive to the shape of the fragmentation functions. 
Our calculations indicate that such measurements could be used to test 
the color-octet mechanism solution to the ip' surplus at the Tevatron. 



1 Introduction 



Recent Tevatron data |T],H on prompt charmonium {ip and production 
at large pt has attracted a lot of attention. 

At large pt, charmonium production is believed to be dominated by par- 
ton fragmentation 0, and the Tevatron data does indeed show a 1/p^ 
behavior consistent with the predictions of a fragmentation process. 

The large c quark mass rric make, however, predictions possible also for 
the magnitude of the cross section. Production of charmonium states can 
be factorized into the hard QCD production of a heavy cc pair within a 
region l/rric, followed by the nonperturbative formation of the charmonium 
state. In , a general factorization scheme was developed to all orders in the 
relative velocity v of the c and c quark in the charmonium state. For S-wave 
states, treated to leading order in f^, this factorization scheme reduces to 
the color singlet model the cc pair must be produced in a color singlet 
state with the same spin and angular momentum as the charmonium state. 
The subsequent nonperturbative formation of the charmonium state can be 
described by a single process independent parameter. Fixing this parameter 
from measurements on the electro-magnetic decay of the ip and ip', predictions 
can be made. 

A color singlet model calculation to LO in of large pt ip' production fall, 
however, more than an order of magnitude below the Tevatron data 
Although this is the largest deviation from LO color singlet model predictions 
seen, it should be noted that large deviations have been observed also at low 

PT Bi- 

A calculation of prompt ip production [^,|^,|T0[ give a similar result, but 
in this case contributions from P-wave charmonium states decaying to ip 
obscure the interpretation of the data. 



2 



Braaten et al. |TT| have proposed that the ip' surplus at the Tevatron could 
be explained by a f suppressed production mechanism, usually referred to 
as the color octet mechanism. The cc pair is produced in a color octet ^Si 
state, and forms the ip' through soft gluon emission or absorption. The new 
nonperturbative parameter giving the "probability" for the color octet cc pair 
to form the ip' is supressed by f ^ as compared to the color singlet mechanism. 
The hard gluon fragmentation production of a color octet cc pair is, on the 
other hand, of lower order in as as compared to the production of a ^5*1 color 
singlet cc pair. 

In addition, the color octet fragmentation process is enhanced by the so 
called trigger bias effect. The fragmentation contribution to the charmonium 
production cross section can be written as a convolution 

dao{p)=y^l dzdai{p/z,^)Di_^oiz,n), (1) 

of parton production cross sections dci and fragmentation functions Di^Q, 
describing the fragmentation of a parton i into a charmonium state O. Since 
the parton production cross sections fall steeply with the parton transverse 
momentum pr/z, fragmentation processes where the charmonium state takes 
a large fraction z of the parton momentum are favored. In a strict LO 
calculation the color octet fragmentation function is proportional to 6{1 —z). 
When leading logs are re-summed, fragmentation functions are softened, but 
the fact remains that in comparison with other production mechanisms the 
color octet mechanism is characterized by large z fragmentation. 

The arguments in favor of the color octet solution are of a qualitative 
nature, given a new nonperturbative parameter one can obviously fit to data. 
It is therefore crucial to find independent means to check the hypotheses by 
Braaten et al.. 

One way would be to look for other processes where the same color octet 
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production mechanism could be important, see eg. Alternatively one 

might look for characteristics of the color octet production mechanism at the 
Tevatron. It has eg. been noted that ip' produced through the color octet 
mechanism are dominantly transversely polarized [pT^ . 

Here we will make use of the fact that a ip' produced in a fragmentation 
process will be part of a jet. We investigate the possibilities to extract more 
information on the fragmentation process by detecting the associated jet. 
At a superficial level one can identify the jet momentum with the parton 
momentum and thus get a measure of the fraction z of the parton momentum 
taken by the ip' . As noted above, large z fragmentation is a characteristic of 
the color octet production mechanism. 

In it was proposed that trigger bias enhanced higher order corrections 
might explain the ip' surplus within the color singlet model. It was shown 
explicitly that trigger bias enhanced higher order corrections to light quark 
fragmentation into rjc dominate over the LO result. For ip and ip' higher order 
calculations are however lacking. If this proposal is correct one would expect 
the fragmentation to be intermediate in hardness between the soft LO color 
singlet mechanisms and the very hard color octet mechanism. 

In order to make a quantitative analysis it is necessary to implement the 
fragmentation processes in a perturbative QCD cascade, allowing the use of 
measurable variables such as 2;jet = PTV^'/PTjet- 

In this paper we present an implementation of the LO color octet and 
singlet production mechanisms relevant for ip' production, but the extension 
to other mechanisms and other quarkonia is trivial. 

In section |^ we describe the implementation of fragmentation production 
of quarkonium states in a leading-log QCD shower generator based on the 
Dipole Cascade Model (DCM), first in general and then more specificly for the 
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LO color octet and singlet production mechanisms relevant for ip' production. 

In section ^ we give some results of our simulations and discuss the possi- 
bilities to distinguish between different fragmentation mechanisms based on 
the shape of fragmentation functions. 

Our conclusions are presented in section ^ which is followed by an ap- 
pendix describing the analytical approximations used to verify our imple- 
mentation. 

2 The Implementation 

The fragmentation production of quarkonia may be viewed as a two step 
process: 1. the hard fragmentation into a color singlet or octet cc-pair, and 
2. the nonperturbative formation of the quarkonium state. 

Here we implement both these steps in conjunction in a hard QCD shower 
generator, by including splitting functions not into cc-pairs but directly into 
quarkonium states. Thus, the nonperturbative formation of the quarkonium 
state is implemented at the perturbative level, before normal hadronization 
sets in. This is motivated by the fact that all hard interactions are already 
included in the splitting functions. The cc-pairs are thus effectively non- 
interacting until the hadronization onset scale. Finally, according to BBL 
factorization the "probability" that a color octet or singlet cc-pair form 
a quarkonium state is given by process independent parameters (matrix el- 
ements). BBL factorization is however fully inclusive and does not describe 
how the cc-pair effects the hadronization of the remaining partons. In the 
color singlet case, the cc-pair is color disconnected from the other partons 
in the event and it is certainly reasonable to let the other partons hadronize 
independently of the cc-pair. The color octet case is more troublesome. In 
our implementation we enforce the radiation of a soft gluon that inherits 
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the color connection of the cc-pair to the other partons. That is, we sim- 
ply implement splitting into the quarkonium state plus a soft gluon. The 
hadronization is done using the stringfragmentation model in Jetset which 
is infrared safe w.r.t. soft gluons, suggesting insensitivity to the exact mod- 
elling of these soft interactions. We will, however, return to this nontrivial 
question in a forthcoming publication. 

In this paper we present an implementation of heavy quarkonium pro- 
duction in the Ariadne program which implements the Dipole Cascade 
Model (DCM) [|TB],[T^ for QCD cascades. 

The DCM differs from conventional parton shower models in that gluon 
emission is treated as dipole radiation from color dipoles between partons. In 
eg. e^e~ annihilation, the first gluon, gi, is emitted from the dipole between 
the initial qq pair. An emission of a second softer gluon can then be described 
in terms of radiation from two independent dipoles, one between the q and 
gi and one between gi and q. Further gluon emissions are given by three 
independent dipoles etc. 

One problem with this model is that only gluon radiation is modelled and 
the process of splitting a gluon into a qq pair is not included in a natural way. 
In |T^ it was shown that this process can be included without leaving the 
dipole picture by taking the ordinary splitting kernel for g qq and dividing 
it in two equal parts describing the contribution from each connecting dipole 
to the splitting separately. This procedure works very well and describes cor- 



rectly the amount of secondary heavy quarks produced at LEP ||T9|] although 
some theoretical questions have been raised ^0 . 



The process of (? O + g will here be treated in the same way as the g 
qq one. The splitting kernel dg^o{z, S), which is related to the fragmentation 
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function through 

Dg^o{z,ix)^ —dg^oiz,S) (2) 

is hence divided into two equal parts assigned to each of the connecting 
dipoles. Assuming a dipole between a quark and a gluon in its center of mass 
system, with the gluon along the positive z-axis, z is then the fraction of the 
total positive light-cone momenta of the dipole carried by the quarkonium. 
The partons in the DCM is always on mass shell, therefore, S is no longer 
the mass square of the split gluon, but is still the mass square of the decay 
products. Four-momentum is conserved by "borrowing" momentum from the 
other parton in the dipole - this is a valid approach in the limit of strongly 
ordered emissions. 

For practical reasons, the generations of quarkonium production is done 
in the variables = + = zS and y = \n , the transverse mass 
and rapidity of the quarkonium in the dipoles center of mass, W is the total 
mass of the emitting dipole. 

For eg. a q-g dipole there are now several competing processes that can 
occur; a gluon can be emitted according to the standard DCM, the initial 
gluon can be split into a qq pair of different flavors and it can be split into 
a quarkonium plus gluons according to several different mechanisms. This 
competition is as usual governed by the Sudakov formfactor mechanism using 
the transverse momentum squared as ordering variable. Hence the probabil- 
ity Pi{k\) of the process i to occur at a scale kj_ is given by 

P,{kl)^d,{kl)e '"'^^'^ (3) 

where di{k'j_) is the splitting kernel with the dependence on rapidity inte- 
grated out. 
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It is now straight forward to implement basically any mechanism for frag- 
mentation production of any heavy quarkonium. In this paper we discuss 
three mechanisms for ip' production: 

• Color-Octet mechanism g — > cci^sf^) soft gluons. A perturba- 
tive transition of a gluon into a collinear cc pair in a color octet ^Si 
state, followed by a nonperturbative transition into ip' involving emis- 
sion and/or absorption of soft gluons. As discussed above we implement 
this process as g tp' + g, with one soft final state gluon . Since the 
gluons are soft, the splitting kernel is essentially a delta-function, 

dg^^>{z,S) = — d{l - z)d{S ~ AmJ. (4) 

In the implementation we give for practical reasons the splitting kernel 
a small width, but we do not in this paper study possible effects of 
the physical width |2l[]. Since the characteristic velocity v of 

the c-quarks in the charmonium state is not that small (f^ ~ 0.2), 
the physical width could effect the event structure and through trigger 
bias enhancement also the normalization of the large pt quarkonium 
production. The Color-Octet matrix element (O|0g ('^ 5*1)10) was taken 
to be 0.0042 GeV^ as given by the fit in Q to CDF data on large pt 
ip' production. 

• Color-Singlet mechanism g — > cci^^Sf^) + 2g ^ ip' + 2g is theoretically 
simpler since the hard splitting color disconnects the cc-pair from the 
rest of the event. The implementation is complicated, however, due to 
the final 3-particle phase space. For simplicity we treat also this process 
as g ^ ip' + g. The error introduced by this simplification, which will 
be investigated in a future publication, may be assumed to be small 
since two gluons nearby in phase space are usually well approximated 
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with one within the framework of Lund string fragmentation. The 
color-singlet gluon fragmentation function is given as a double integral 
in 



2^ . Performing one of the integrals we find an analytical expression 
for the fictitious g '4'' +9 splitting kernel. The color singlet parameter 
|/2(0)p relevant also for the last production mechanisms is taken to be 
0.231 GeV^. 

Color-Singlet mechanism for c ip'+c. Here the whole splitting kernel, 



given in [Q, is assigned to the only dipole connected to the c quark. 
The scale is the running as is, as in the other cases, taken to be m^. 

3 What comes out 



The Ariadne program is interfaced to Pythia and Jetset gS] for 
handling of the hard sub-process and hadronization respectively. We use 
this to generate ip' production events at the Tevatron. In fig. [1| we see the 
resulting pT^p' spectrum for the different production mechanisms. Using the 
same parameters as in [^, we reproduce approximately their result where 
the Octet mechanism totally dominates the cross section for high pTip', and 
we also reproduce the measured Pn^'-^^i+^i- spectrum 

The difference between the Octet and Singlet mechanisms is larger in our 
case than in even after taking into account the small variations in the 



parameters used. The generated spectra from Ariadne does, however, agree 
with the analytical approximation in the appendix at large pt- At small pt 
we do not expect these to agree because there the simple assumption in the 
analytical calculation that the gluon spectrum behaving like breaks 
down. Also effects of initial-state QCD radiation, implemented through the 
so-called Soft Radiation Model p6[ in Ariadne, are absent in the analytical 
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Figure 1: The transverse momentum distribution of ip' ^~^fi~ production 
at the Tevatron for different production mechanisms compared to CDF data 
Full line is the color octet mechanism, dotted line is the sum of the 
color singlet mechanisms, and the long-dashed and short-dashed lines are the 
gluon and c- quark initiated singlet mechanisms respectively. The histograms 
were generated with Ariadne and the smooth lines are from the analytical 
approximation in the Appendix. 



calculation. In addition, Ariadne implements a more complete treatment 
of the parton evolution, in contrast to the analytical calculations, which, like 
earlier numerical calculations, uses the approximative though asymptotically 
correct AP evolution. The full evolution equation and a discussion of the 
approximation made in the use of the AP evolution equation can be found 
ini. 

In fig. 1^ we show the generated 2;-spectrum for 20 < pti(>' < 25 GeV, 
where z is now PtxP'/pt and pT is the transverse momentum in the hard 
subprocess. It is clear that the z-spectrum is much harder for the octet 
mechanism as expected. This is also seen in fig. |^, where the average z is 
shown as a function of pr^p'- Here we also find a good agreement with the 
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Figure 2: The z- distribution for 20 < ptii,' < 25 GeV as generated by Ari- 
adne for different production mechanisms. 
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Figure 3: The average z as a function of pT^p' for different production mech- 
anisms. The smooth lines are from the analytical approximation in the ap- 
pendix and the (due to low statistics) broken lines were generated with Ari- 
adne. 
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Figure 4: The z distribution for 20 < pTip' < 25 GeV as generated by Ari- 
adne using the color octet mechanism. Full line is the generated spectrum, 
long-dashed line is after fragmentation using a cone radius of 1 to define the 
associated jet, and short-dashed is the same with the addition of the under- 
lying event according to the multiple scattering model in Pythia. 



analytical approximation at large Ptv^'- 

With the availability of ip' production in an event generator we can now 
also study hadronization effects. To study eg. the 2;-spectrum experimentally 
one could sum up the transverse energy -Erjet in a cone in the pseudorapidity- 
azimuth plane around the detected ip\ defining z-^et = PT^p'/ETjet- In fig. ||we 
compare the 2;-spectrum for such a procedure with the generated spectrum 
in fig. 1^ for the octet mechanism. We note that just adding hadronization, a 
smearing of the jet is introduced and some of the jet energy falls outside the 
cone, which makes the spectrum considerably harder. Adding also multiple 
interactions to model the underlying event ||2^ in Pythia makes it softer 
again. In fig. |^ we see that the differences in hardness between the octet 
and singlet mechanisms remains after all hadronization effects although the 
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Figure 5: The z distribution for 20 < pT^p' < 25 GeV as generated with 
Ariadne with different production mechanisms, after hadronization and in- 
clusion of multiple interactions. 

shapes are somewhat distorted. 



We have implemented different quarkonium production mechanisms in the 
Ariadne event generator. Besides giving a more accurate treatment of the 
initial- and final-state parton evolution than recent analytical approxima- 
tions, such a generator is an essential tool for understanding experimental 
corrections and hadronization effects in case one wants to study non-inclusive 
aspects of high p^ onium production at e.g. the Tevatron. 

In interpreting the results it is important to keep in mind the status of 
the implementations of the different production mechanisms: 

• The implementation of the color singlet c ^ ip' + c fragmentation is 
unproblematic and trustworthy. 



4 Conclusions 
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• The color singlet g ^ ijj' + gg fragmentation has been implemented 
ets g ^ ip' + g. The error introduced by this simplification is believed 
to be small, and in the future the full g ^ ip' + gg splitting will be 
implemented. 

• The color octet g — > cc{^Si^^) — > soft gluons, has been implemented 
as g ^ '4>'-\- one soft gluon, and the effects of the physical 0{v'^) width 
of the fragmentation function has not been investigated. In the future 
the sensitivity of the results to the modeUing of the soft interactions 
will be studied. 

The production mechanisms we have implemented here are, after a sim- 
ple change of normalization relevant also for production. In this case it 
is necessary, however, to take into account also contributions from P-wave 
states decaying to -0. It is therefore important to implement also P-wave 
production mechanisms. 

There are also other experiments where the fragmentation production 
of quarkonium states is important. One interesting example is decays at 
LEP. Here the situation is quite different as compared to large pt production. 
There are more quarks than gluons available for fragmentation, and there is 
no trigger bias enhancement. Consequently other production mechanisms 
might be dominant and measurements would be more sensitive to the low z 
shape of fragmentation functions. 

We have shown that it might be possible to measure the ^-distribution 
of the fragmentation function at the Tevatron and from that learn about the 
production mechanism. 

This result is of course dependent on our very crude approximation of the 
experimental environment. But with the implementation presented here, the 
determination of optimal choices of jet definitions and experimental cuts is 
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clearly facilitated and it does not seem unlikely that the z spectrum, if mea- 
sured at the Tevatron, could be used to better understand the quarkonium 
production mechanism. 

Acknowledgement. We are greatful for discussions with Paul Hoyer, 
Michelangelo Mangano and Mikko Vanttinen. 
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Appendix 



Here we present a very simple calculation that exhibits the gross features of 
charmonium production at large Pt- 

We start out from a rapidity integrated version of eq. (|l]) 

We choose the factorization scale as fi = pr/z, and note that large logs 
log(4m^//i^) can be resummed by AP evolution ||28| of the fragmentation 
functions. We calculate the parton production cross sections as convolutions 
of hard subprocess cross sections (oc p^^) and parton distributions, and find 
that the parton production cross sections fall roughly as p^^. Consequently 
we take dai/dpij^ ~ Ci/pij^^ in eq. @ and find 



^ - T.£dz^^'D.-^A^,p^T/z)^T.^DrUpr). (6) 
apT i JO Pt i Pt 



Here Df^^,{jj,) = Jq dzz^Di^^'{z, jj) is the sixth moment of the fragmentation 
function. In eq. (^) we have neglected the evolution of the fragmentation 
functions between fi = pt and Pt/ implicitly assuming that the dominant 
contributions to the integral in eq. come from large values of z. 

We find ^ 7 ■ lO^GeVSb, and ^ 4 ■ lO^GeVSb, to LO in a,. The 
gluon production cross section is two orders of magnitude larger than the 
quark production cross section. Within the color singlet model to LO in 
«<(, the c quark fragmentation function is, however, much larger than the 
gluon fragmentation function. We keep therefore both gluon and c quark 
fragmentation in eq. (P). 

In Q mixing effects were taken into account by including NLO corrections 
to the hard scattering charm production cross section. Here as well as in our 
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QCD shower generator implementation, we include mixing explicitly and use 
the LO charm production cross section to avoid double counting. 

For moments (n > 1) the AP evolution equations can be solved analyti- 
cally. The solution has the form: 

-,(") 



4"' 



+ 



r)(") 



+ 



with p(")5(") = e("^(l — e^"^). For the moments n relevant here, the mixing 
parameters p*^"-', 5*^"^ and e*-"-* are all small, and we may put p*^"-* and e*^"^ 
to zero. Within the color singlet model the terms in eq. (0) are 

important however since as noted above the c quark fragmentation function 
is much larger than the gluon fragmentation function. For simplicity we 
choose the initial scales for charm and gluon fragmentation and the mixing 
onset scale all to 2mc and ignore flavor thresholds (refinements are simple). 
Combining eq. (|^),(^ and (||) we find 



I- 



(6) 

Dfi^,(2m,)-5(«)Z}f^,(2m,)) + 

(6) 

(6) 



5 ^, (O^ \ c^'/' 



\as{2m. 



D':L,{2m^). (9) 



According to the hypotheses by Braaten et al. []TT]| the color octet mechanism 



gluon fragmentation completely dominates ip' production. We may then of 
course ignore charm fragmentation all together in eq. (^. 
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Similarly we can calculate the mean value of z in the convolution in 
eq. (|^). Superficially identifying z with PT^i/pT^et, this gives us a hint of 
the possibilities to use the form of fragmentation functions to distinguish 
between different production mechanisms. More importantly, it can be used 
to check our QCD cascade implementation of quarkonium production. 

For convenience we give numerical values for the moment evolution pa- 
rameters '^i"^ 4 (5) quark flavors, as well as the values we have 
used for the moments of the color singlet fragmentation functions D^^^i and 



n 


5 


6 


7 


Sin) 


0.030 (0.028) 


0.023 (0.022) 


0.018 (0.017) 




2.18 (2.46) 


2.46 (2.76) 


2.68 (3.00) 


4") 


0.97 (1.05) 


1.08 (1.17) 


1.17 (1.27) 




1.2 • 10-5 


9.7-10-^ 


7.7-10-^ 




1.8-10-^ 


1.3 ■ 10-^ 


1.1 ■ 10-7 



For the moments of the color octet gluon fragmentation function we have used 
the value D^^}^^,{2mc) ~ 4.2 • 10"^, independent of the moment n. 
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